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New data from WMAP have appeared, related to both the fractional energy density in relativistic 
species at decoupling and also the primordial helium abundance, at the same time as other inde- 
pendent observational estimates suggest a higher value of the latter than previously estimated. All 
the data are consistent with the possibility that the effective number of relativistic species in the 
radiation gas at the time of Big Bang Nucleosynthesis may exceed the value of 3, as expected from 
a CP-symmetric population of the known neutrino species. Here we explore the possibility that new 
neutrino physics accounts for such an excess. We explore different realizations, including neutrino 
asymmetry and new neutrino species, as well as their combination, and describe how existing con- 
straints on neutrino physics would need to be relaxed as a result of the new data, as well as possible 
experimental tests of these possibilities. 



INTRODUCTION 

One of the remarkable successes of Big Bang Nu- 
cleosynthesis (BBN) is the correct prediction of the 
overall magnitude of the light clement abundances as 
a function of two fundamental parameters: the baryon 
density in the universe, and the number of light neu- 
trino species. Until recently BBN provided the only 
direct handles on these two fundamental parameters 
(i.e. [THH). However with the discovery of primor- 
dial anisotropics in the Cosmic Microwave Background 
(CMB) radiation, the CMB has become a remarkable 
precision laboratory to constrain fundamental param- 
eters in particle physics and cosmology. It can now be 
used to test various ideas associated with BBN, and 
consistency checks can be applied to probe new physics. 

It is therefore of some interest that recent results 
from both of these areas suggest the possibility the 
some new physics beyond the standard model may be 
at play. Izotov and Thuan 4] have recently analyzed 
the primordial helium abundance. For the primor- 
dial 4 He mass fraction, Y p , they find Y p = 0.2565 ± 
0.0010(stat.) ±0.0050(syst.), which is higher, at the 
2a level, than previous measurements (see e.g. [5] and 
references therein.) Note that this value is in good 
agreement with another recent estimate [6] although 
their quoted error is far smaller. At the same time 
the new WMAP 7 year analysis [7] suggests both a 
high value of Y p (with larger error bars), and inde- 
pendently a somewhat high value of the N e f f , the ef- 
fective number of relativistic neutrino species present 
during last scattering (N eff = (p Tel - p 7 )/p Utbmia , 
where p„ tbmm = (77r7l20)(4/ll) 4 / 3 T 4 ). Specifically, 
the WMAP 7 measurement is N eff = 4.34 + 0.86-0.88 
[Sj, about 1.4 o higher than the standard contribution 
of the known neutrino species, N e ff= 3. 

While the uncertainties in these estimates remain 



large, if confirmed experimentally in the future, a sce- 
nario with high Y p and high N e f / would be an interest- 
ing challenge for particle physics. While multiple, un- 
correlated, phenomena could explain it, both the Izo- 
tov and Thuan and the WMAP 7 results are clearly 
consistent with a single physical cause: that number of 
helicity degrees of freedom in the radiation gas at the 
time of BBN might have been significantly higher than 
3. The simplest and most conservative possibility for 
this involves new neutrino physics. 

Within the context of neutrino physics alone, two 
different possibilities arise. One can consider mecha- 
nisms which affect both the energy density of neutrinos 
during BBN and the weak interaction rates that deter- 
mine equilibrium nuclear abundances, such as would 
be the case for a neutrino-antincutrino asymmetry. As 
a more minimal option, a correlation between two phe- 
nomena that are very separated in cosmic time-an el- 
evated helium abundance after BBN (T ~ 0.2 MeV) 
and extra energy density in relativistic degrees of free- 
dom at matter- radiation decoupling (T ~ 0.1 eV)- 
could have a common origin in an overabundance of 
weakly interacting and light (sub-eV) mass particles at 
the time of BBN. During BBN this overabundance en- 
hances the energy density, which, as we have alluded, in 
turn increases the expansion rate of the universe. This 
causes the weak reactions to freeze out earlier, result- 
ing in a higher neutron-to-proton ratio and therefore 
a higher Y p [H I5HT5] . Depending on their coupling 
to matter, the same particles could also affect Y p by 
contributing to the weak reaction rates. At matter- 
radiation decoupling, they could still be relativistic and 
therefore contribute to N e ff as measured by WMAP7. 
Other possibilities will also be briefly discussed. Previ- 
ous strong constraints on both weak interaction physics 
and neutrino flavors need to be relaxed as a result of 
the current data, allowing possibilities that had previ- 
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FIG. 1: Y p contours in the £„ c and N e ff parameter space 
assuming neutrino flavor equilibration (£„ e = = £,v T )- 
The horizontal light (yellow) band corresponds to the la 
WMAP 7 year result. The black contours show a range of 
calculated values of Y p given model independent inputs of 
£„ e and N e s. The shaded (colored) vertical bands mark the 
Izotov and Thuan la, 2a, and 3a ranges of Y p . The bottom 
black curve shows the contribution to N e ff from neutrino 
asymmetries alone. 

ously been considered as ruled out. 

Clearly, suitable candidates should fit a number of 
conditions: (1) they should couple to matter and radi- 
ation strongly enough to be produced by thermal pro- 
cesses before BBN but (2) should not contribute too 
many extra degrees of freedom to the radiation gas, 
constrained in turn by measurements of Y p . Further- 
more, (3) they are constrained by the existing cos- 
mological bounds on the density of light extra de- 
grees of freedom coming from a combination of data 
from the CMB, Large Scale Structure (LSS), Lyman 
Alpha Forest, and Baryon Acoustic Oscillations [5]: 
Q„h 2 < 0.006 (95% CL). We consider specific scenarios 
and constraints in the following sections. 

NEUTRINO ASYMMETRIES AND DECAYS 

An overabundance of neutrinos with respect to anti- 
neutrinos or vice-versa, L v = (n u — n p )/n 7 , is de- 
fined by a non-zero degeneracy parameter, £ : L v = 
ir 2 /(12((3))(T„/T) 3 (£ + C 3 A 2 )- The total change in 
the effective number of relativistic species resulting 
from asymmetries in each flavor, £„ Q , is given by 




In most theoretical scenarios, lepton and baryon asym- 
metries are enforced to be of the same order by 
sphalerons [TS], so that L v ~ 10~ 10 — 10~ 9 . However, 
several scenarios have been proposed in which a large 
lepton asymmetry can be generated while preserving 
a small baryon asymmetry, using e.g., GUT models, 
the Affleck- Dine mechanism, Q-balls, resonant oscil- 
lations, etc. [TTHUj. Therefore, here we assume L v as 
independent from the baryon asymmetry and consider 
only direct constraints on it from neutrino physics. 

While asymmetries in all flavors contribute to an 
increase in energy density, only an asymmetry in the 
electron flavor influences the weak neutron-proton in- 
terconversion processes. For this reason, the sensitivity 
of BBN to £ u<j is remarkably high: |£„J <J few 10~ 2 
is needed for compatibility with measured abundances 
(see e.g. [2 [H 113 [23H22])- This applies also to the 
asymmetries in the other flavors at the time BBN, since 
oscillations should produce an at least approximate fla- 
vor equilibration before BBN [2"8H31| . 

Under such strong constraint, neutrino asymmetries 
alone generally cannot account for a AN e ff ~ 1. An 
interesting exception is the somewhat fine-tuned sce- 
nario of initial (pre-equilibration) flavor asymmetries 
that are large and opposite in sign. After equilibra- 
tion, a surviving AN c g ~ 1 can be realized, together 
with sufficiently small asymmetries that satisfy BBN 
bounds [35]. This reopens the possibility of having, 
at BBN, virtually any combination of £„ a and energy 
density. In general, asymmetries could coexist with 
other effects (e.g., a sterile neutrino, see next section) 
that could independently increase N e ff. Therefore an 
analysis that treats asymmetries and energy density 
as independent is necessary to find the most general 
constraints on both. 

Here we perform such a study, using a modified ver- 
sion of the Kawano/ Wagoner BBN code described in 
detail in Ref. [T3] [33]. In Fig. 1 we illustrate the 
interplay between asymmetries and N e ff by plotting 
the Y p abundance yield isocontours in the N e ff - £„ c 
parameter space. This figure shows calculations for 
model-independent inputs of N e ff over a wide range 
of neutrino asymmetries, where we have adopted the 
condition of neutrino equilibration of asymmetries, so 
that the allowed range of asymmetries is small, and the 
direct effect of such asymmetries on N e ff is minimal 
(as displayed in the lower curve, which shows the extra 
direct contribution to N e ff from such asymmetries). 
The horizontal band for N e ff corresponds to the la 
WMAP 7 year result quoted earlier. 

The BBN code used to make Fig. 1 differs from oth- 
ers in a number of ways, mostly in the treatment of 
the weak processes. It allows for calculations that in- 
clude both the effects from higher relativistic degrees 
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of freedom and neutrino asymmetries, which is not the 
case for the Kawano/Wagoner code. This code gives 
a standard BBN Y p yield that is ~ 0.004 lower than 
other calculations due to the full numerical integra- 
tion of each weak reaction rate. This figure is not 
intended to provide new constraints on, or a best fit 
for, neutrino asymmetries and/or N e ff. It is simply 
a model-independent tool to illustrate the fact that a 
wider range of allowed N e ff parameter space loosens 
the BBN limit on the neutrino asymmetries or lepton 
numbers. 

This figure demonstrates explicitly the (known) fact 
that the neutrino electron degeneracy parameter has 
significant leverage on Y p . Positive £„ drives the neu- 
tron destruction process, v e + n ^ p + e~ forward 
and and Pauli-blocks the reverse neutron production 
reaction resulting in a lower Y p . Negative £„, corre- 
sponding to an overabundance of anti-neutrinos, drives 
anti-neutrino capture forward v e + p # n + e + and 
suppresses the reverse proton production reaction re- 
sulting in more helium. This figure displays explicitly 
the novel recognition that to obtain a given isocon- 
tour for Y p , one can increase while also increasing 
N e ff, and it provides a quantitative estimate of the 
interplay between these effects. This interplay has not 
been directly examined quantitatively before, although 
a correlation between the two parameters was noted in 
a fit to BBN data [55] and in an analysis of the baryon- 
to-photon ratio[14). 

The figure also shows how the constraints on ^ 
change with a change of N e f f and can thus be relaxed 
compared to previous limits. For example, if we al- 
low N e ff in the la WMAP7 interval and Y p in the 3a 
interval of Izotov and Thuan we get an approximate 
allowed range of 

- 0.14 <> £ 0.12, (2) 

larger than the usually quoted constraint —0.04 < 
£v e < 0.07 [53]. While these results present a model- 
independent exploration of parameter space, with each 
point involving a full BBN code calculation of Y p for a 
given value of N e ff and £, whether any specific point 
in parameter space is actually realizable however, will 
depend upon specific model building issues. 

We conclude this section by briefly mentioning an- 
other possibility for adding extra relativistic energy 
density both during BBN times and during the matter- 
radiation equality epoch: extra particles that are un- 
stable. To produce higher Y p , these particles must con- 
tribute to the relativistic energy density during the 
"weak freeze out" period in BBN, implying masses 
rn < 1 MeV. They would then be required to decay 
by the time of matter-radiation equality (T ~ 1 eV), 
where the CMB measurements infer extra relativistic 



degrees of freedom (see for example [541 I35] l. 

However, additional cosmological constraints imply 
that the physics of such particles is so finely tuned as 
to be implausible. Their decay can only be into neutri- 
nos so as not to produce high-energy photons which 
result in subsequent deuterium photo-disassociation 
[55] . Furthermore, the additional particles must de- 
cay quickly enough so that they don't subsequently 
dominate the energy density of the universe once the 
temperature falls below their mass. The dual require- 
ments of being primarily weakly interacting and also 
decaying within the appropriate time window are ex- 
tremely difficult to satisfy. 

STERILE AND RIGHT HANDED NEUTRINOS 

The minimal scenario to explain both high N e f f and 
high Y p with neutrino physics is a light sterile neutrino 
(see also [37 ). Indeed, sterile neutrinos easily fit the 
conditions we have outlined above: by definition they 
are weakly interacting, they can be produced before 
BBN, and they are allowed, by laboratory and astro- 
physical bounds, to be of sub-eV mass. 

Specifically, if the neutrino masses come from a 
See-Saw-like mechanism, active-sterile oscillations arise 
naturally due to the mixing of active neutrinos with 
the charge-conjugate of one or more right handed neu- 
trinos. Depending on the mixing and masses, the in- 
terplay of oscillations and collisions can populate the 
sterile neutrinos before BBN (see e.g., [35] and refer- 
ences therein). While in minimal See-Saw models the 
sterile neutrinos are too heavy to contribute to N e ff, 
several non-minimal scenarios (e.g., [39[ I40j ) include 
sterile neutrinos lighter than an eV. 

Detailed studies exist on the effect of one ster- 
ile neutrino (v s from here on) on Y p and N e ff (e.g. 
[5TI [HJ [32] ) . Here we consider one such case, assuming 
a hierarchical spectrum with the predominantly sterile 
state being the most massive, Am^ 7$> \Am%i\. Fig. [2] 
(adapted from the results of [32]) refers to the specific 
case in which v s mixes with 2/„ and v e , as needed to in- 
terpret the LSND anomaly [33] . It is however represen- 
tative of the general situation. The four-neutrino mix- 
ing scheme can be reduced to an effective two-neutrino 
mixing with parameters Am LSND ~ Arn^j = m\ — m\ 

and Qlsnd — S es 9 ps [31] (Figure [2] refers to the case 

1/2 .... 
® lsnd — ^es — Ops, which is conservative in that it 

generally corresponds to minimizing the high Y p re- 
gion, see [4"2"].l. 

We consider both sterile neutrinos with and without 
lepton asymmetry, and we comment on the case of zero 
asymmetry first. In the plot we highlight the region of 
interest: the area (light shaded in the plot) where v s 
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FIG. 2: Limits on the sterile neutrino oscillation parameters 
Am 2 LSND ( m eV 2 ) and 8lsnd from several neutrino experi- 
ments and from the cosmological bound Q, v h 2 < 0.006. The 
shaded regions are allowed. For zero asymmetry, the dark 
(light) shaded area corresponds to Y p < 0.258 (Y p > 0.258). 
The dashed line represents the boundary of the Y v > 0.258 
region for asymmetry L v = — 10~ 6 . The region favored by 
the LSND data, now almost entirely excluded, is shown as 
well. See text for details. 

is produced abundantly prior to BBN, thus causing a 
high Y p via its contribution to the energy density, and 
contributing to N e ff at matter radiation-decoupling. 
Specifically, the region corresponds to Y p > 0.258 and 
N e ff= 3.8 — 4, i.e., a nearly or completely populated 
sterile state. It is bounded from below by the "ther- 
malization line" , where the v s production rate is com- 
parable to the cosmic expansion rate [H]. and spans 
more than an order of magnitude in each parameter, 
extending down to sin 2 20lsnd ~ 10 -5 . The region 
is constrained in mixing by several terrestrial experi- 
ments (mainly Karmen, Bugey, SuperK, CDHS 0SH1H] 
[55] and in Am^ SIVD by the cosmological bound on Q„. 
For N e ff= 4, this bound gives: 

Amli & 0.28 eV 2 m 4 < 0.53 eV , (3) 

assuming, conservatively, mi ~ m2 <C ~ 0.05 eV, 
as given by oscillation data for the normal mass hierar- 
chy [5] . We stress that the area we are considering was 
interpreted as excluded by BBN until recently. Our 
new perspective reopens this possibility. 

Results similar to those in fig. [2] are obtained for 
other active-sterile mixing scenarios, such as those in 
which v s mixes with one active flavor or with one neu- 
trino mass eigenstate. The mixing of v s with is the 
least constrained because of the strong constraint on 
the v e component of v% |41j . 

In the presence of lepton asymmetry, the produc- 
tion of sterile neutrinos via oscillations is suppressed 



[HI H5H5T] . This effect is due to the term in the po- 
tential describing neutrino-neutrino forward scattering, 
which suppresses the active-sterile mixing for neutri- 
nos or antincutrinos and is zero for a symmetric neu- 
trino population. Fig. [2] shows how the high Y p region 
changes with the increase of L v which is assumed to 
be equal for all flavors. As a result of the suppression, 
with lepton asymmetry the region of high Y p is reduced 
to a smaller area at high Am 2 LSND (fig. [2]) and eventu- 
ally disappears for L„ ~ 10 -5 [42] as the mass required 
to populate the sterile neutrino becomes excluded by 
the bound on £l v . Thus, a single light sterile neutrino 
can give only a subset of the region in our fig. [I] This 
subset has 3 <N e ff < 4, with N e ff~ 4 being realized 
only for £„ e ~ 0. 

The conclusions on the suppression of the production 
of v s are generally true for a wide range of lepton asym- 
metry, L v ~ 10~ 5 — 1, provided that the asymmetry is 
constant over the characteristic time scale of the sterile 
neutrino production. For L v ~ 0.1 — 1 the presence of 
the sterile state can modify Y p through a modification 
of the spectra of the active states, while not affecting 
N e ff [2U[S5J[S2]. Therefore we do not consider this sce- 
nario here. A more diverse phenomenology is expected 
if the asymmetry varies over the time of v s genera- 
tion: active-sterile oscillations can actually generate an 
asymmetry in the active flavors that can survive and 
affect the weak reaction rates [M] [55] ■ Although an 
updated analysis on this is not available, from existing 
studies |51j we infer that if the sterile neutrino con- 
tributes substantially to N e ff, and therefore its mass 
is below the LSS bound, the generated asymmetry is 
L v ^ 10~ 2 , not sufficient to impact Y p via BBN re- 
action rates. Therefore, this effectively reduces to the 
case L v = 0. 

If the neutrino mass arises from a Dirac mass term 
only, the right handed neutrino(s) associated to it are 
not produced via oscillations from the active states and 
therefore do not play the role of sterile neutrinos as de- 
scribed above. Still, a number of models exist in which 
right handed neutrinos are of sub-eV mass, and cou- 
ple to the Standard Model particles strongly enough to 
be populated substantially prior to BBN [TlJ [56j [57] . 
A detailed analysis in the context of an E 6 symmetry 
|57j shows how, indeed, increased N e ff and Y p are ex- 
pected due to the right handed neutrino production. 
Compatibility with BBN translates into lower limits 
on the mass of the Z' . These limits will be relaxed for 
increased N e ff and Y p \ the new expected mass range 
will be close to or overlap with the limits from SN1987A 
[55] . and therefore measurements from a future galac- 
tic supernova could test such class of models, as we 
describe below. 
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EXPERIMENTAL SIGNATURES 

A sterile neutrino with parameters in our high Y p re- 
gion has a number of implications for future detectors. 

Beta decay and neutrino-less double beta decay ex- 
periments, designed to measure the neutrino mass, 
would probe a fourth light mass state that mixes with 
the electron neutrino (see e.g. [59 ). The next gener- 
ation of reactor neutrino experiments will probe this 
region beyond the existing limits, to an extent that 
depends on the specific model and on #13 [60]. Neu- 
trino beams will also allow to search for sterile states, 
probing different parameters depending on their en- 
ergy and baseline 61j. For a ~ 10 GeV beam and 
Amfj <, 0.1 eV 2 , a baseline of L ~ 2irE/Aml 1 £ 100 
Km is required. Signatures of a sterile state would 
be disappearance of the active flavors and anomalous 
differences in the oscillation pattern of neutrinos and 
antineutrinos due to refraction in the Earth. A recent 
example of the latter involves a sterile neutrino with 
parameters in our region of interest (high Y p and high 
N e ff) [22], and is still interesting in the light of the hint 
of neutrino-antineutrino differences at MINOS [55] . 

It was observed [42] that the suppression produced 
by an asymmetry relaxes the cosmological bound on 
lS.m 2 LSND , thus allowing parameters that explain the 
LSND anomaly. As fig. [2] clarifies, however, these pa- 
rameters correspond to low Y p and N e ff~ 3 (the ster- 
ile state is not populated). Therefore, if MiniBOONE 
confirms LSND and high N e tf and Y p are established, 
less minimal scenarios, beyond a CP-symmetric system 
of four neutrinos, would have to be considered. 

Important astrophysical tests of a light sterile neu- 
trino would come from atmospheric neutrinos and a fu- 
ture supernova neutrino detection |69j Higher precision 
measurements of atmospheric neutrinos would extend 
the currently probed region of parameters. The detec- 
tion of 0.1-1 TeV atmospheric neutrinos at IceCUBE 
would facilitate searches of sterile neutrinos in the 
higher Am 2 range of our region, for which the active- 
sterile mixing is enhanced by matter effects [Ml |SS] . 

In a supernova, a sterile neutrino in our high Y p re- 
gion could be produced via resonant oscillations and 
cause a suppression of the active neutrino signal. Due 
to partial violation of adiabaticity [3"51 |4"T] . the suppres- 
sion would be moderate, at the level of tens of per cent. 
It will be detectable with a future galactic supernova 
if precise theoretical predictions of supernova neutrino 
emission are available. 

Compared to a sterile neutrino, a right handed neu- 
trino would be more difficult to test experimentally, 
because it requires non-oscillation tests. A possibility 
is to look for the new gauge bosons that couple to the 
right handed states, at man made or cosmic acceler- 



ators. Among the latter, the neutrino burst from a 
galactic supernova would be sensitive to the the pro- 
duction of right handed states via their effect on the 
cooling rate, as mentioned previously. Constraints on a 
right handed neutrino will be highly model-dependent. 

If a high N e f t and high Y p are confirmed and other 
data exclude extra neutrino species that can be popu- 
lated before BBN, one would have to consider more 
finetuned scenarios like large and opposite neutrino 
asymmetries in the different flavors. This would indi- 
cate CP violation in neutrinos, with profound implica- 
tions on our understanding of the lepton sector. While 
to uniquely trace back to a model that generates oppo- 
site asymmetries might not be possible, a large number 
of mechanisms that predict more natural, comparable, 
asymmetries would be disfavored. Note, however, that 
it would be very difficult to confirm this scenario inde- 
pendently, because direct CP violation experiments in 
the neutrino sector would not probe the CP violating 
effects of relevance during BBN. Therefore precision 
cosmology alone might be required to more firmly con- 
strain the neutrino sector. 

Precision cosmological tests may be possible in the 
not-too-distant future. Direct kinematic constraints 
on neutrino masses from large scale structure tests, 
in particular from probes of galaxy clustering via the 
Sunyaev-Zeldovich effect on CMB, as currently being 
explored in the South Pole telescope [Bol IST] could im- 
prove existing mass constraints by a factor of 2-5. At 
the same time, the Planck satellite will improve con- 
straints on N e f f by increased sensitivity to high-/ CMB 
anisotropics. Finally, direct measurements of primor- 
dial Y p might be possible in the more distant future 
if spectral sensitivity to the CMB can be improved by 
several orders of magnitude. 

To conclude, if the tentative new results on N e ff 
and Y p are confirmed, they open up new possibilities 
for neutrino physics which may be accessible in the 
future by ground-based experiments and astrophysical 
probes. Specifically we find: 

• A new relaxed constraint on possible electron 
neutrino asymmetry —0.14 ^ £„ e ^ 0.12, 

• if N e f f m 4, a bound Am 2 x < 0.28 eV 2 , m 4 < 
0.53 eV on sterile neutrino masses under fairly 
conservative assumptions about a mass hierar- 
chy 

• a new quantitative relation between the effects of 
possible independent variations in £ Ue and N e ff 
on Y p . 

• a new set of possible astrophysical and experi- 
mental signatures that might further probe these 
scenarios 
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